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^ . Abstract 



^ We give a perturbation theory of neutrino oscillation in the Earth. The perturbation 
g_) ' theory is valid for neutrinos with energy E ^ 0.5 GeV. It is formulated using trajectory 
r^, ' dependent average potential. Non-adiabatic contributions are included as the first order 
effects in the perturbation theory. We analyze neutrino oscillation with standard matter 
effect and with non-standard matter effect. In a three flavor analysis we show that the 
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CN I perturbation theory gives a precise description of neutrino conversion in the Earth. Effect 



of the Earth matter is substantially simplified in this formulation. 
PACS: 14.60.Pq, 13.15. +g 



^ ■ 1 Introduction 

One of the major concerns in neutrino oscillation experiments is the effect of Earth matter 
^ ' in neutrino flavor conversion[ll |2]. Because of the complicated matter density proflle in 
■ the Earth it is a challenge to flnd a simple formula which can describe precisely neutrino 
conversion in the Earth. We rely a lot on numerical computation which does not offer 
us enough insights. Few years ago, an analytic and precise formulation was obtained for 
solar neutrinos, i.e. for E ^30 MeV [SI II]- It is the purpose of the present paper to flnd 
a good formulation for higher energy neutrinos, i.e. for E ^ 0.5 GeV. Previous works on 
neutrinos of this energy range include O El [TJ [U [9l [TOl [HI [121 [ISl • 

It will be shown that the approach using trajectory dependent average potential is 
very attractive. This approach was used in [13] in discussing solar neutrinos. It was used 
in [15] in a 2z/ analysis for higher energy neutrinos {E ^ 0.5 GeV). In the present paper 
we present a perturbation theory for three flavors of neutrinos. The perturbation theory 
uses trajectory dependent average potential. We do perturbative expansion using the 
deviation of potential around the average. We analyze the perturbation theory using the 
density proflle of the Preliminary Earth Model(PREM) We study the perturbation 
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theory for the case of standard matter effect and the case of non-standard matter effect 
with non-standard Neutral Current (NC) interactions. We show that the perturbation 
theory works very well for both cases. The theory is valid for neutrinos with ^ 0.5 
GeV. 

The perturbation theory presented is very useful for long baseline neutrino experiments 
[m im [m [201 121]. it says for a fixed baseline ( < 6000 km) the standard matter effect 
is a one parameter fit. Non-standard matter effect is also greatly simplified. The pertur- 
bation theory is of general interests to other types of neutrino sources, e.g. for studying 
atmospheric neutrinos, cosmic neutrinos from the galactic or extra-galactic sources. The 
present paper is organized as follows. In section [2] we present the perturbation theory of 
three flavor of neutrinos. In section [3] we study neutrino oscillation in the case of standard 
matter effect. In section H] we extend the discussion on neutrino oscillation to the case of 
non-standard matter effect with non-standard neutrino NC interaction. We discuss and 
summarize in section O 

2 Formulation of f oscillation in the Earth 

We consider oscillation of three flavors of neutrinos: = (z/g, t'^, z/,-). The evolution 
equation is 

i£^{x)=H{x)tP{x), H{x) = Ho + V{x) (1) 

Ho = diag{0. Ami, A^li) U\ (2) 

where V{x), a 3 x 3 matrix, is the potential term accounting for the matter effect. V{x) 
takes different form if neutrino interaction is different. U is the 3x3 neutrino mixing 
matrix in vacuum. U is parameterized using standard parameters 612, 61^, 623 and 613, 
the CP violating phase. 

In solving the problem of neutrino oscillation we introduce the average potential term 

V: 

V=]r f dx V{x), (3) 
L Jo 

where L is the length of neutrino trajectory in the Earth. Note that V depends on the 
trajectory of neutrino in the Earth. It is averaged over potentials along the trajectory 
and is not averaged over potentials of all points in the Earth. Using V we introduce H 
and the mixing matrix Um in matter which diagonalizes H: 

H = Ho + V, (4) 
HU^ = U^^A, A = diag{Ai,A2,A3}. (5) 
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2^A*(i = 1, 2, 3) are three eigenvalues of H. < < A^ is satisfied in our convention. 
Um is parameterized using parameters 6*^, 6*]^^, 6*^ and the CP violating phase 5]^. Note 
that we work in the convention of Hamiltonian introduced in and (jl]). In typical 
cases where 3x3 evolution problem can be reduced to 2 x 2 evolution problem, e.g. when 
hierarchy in eigenvalues are present, eigenvalues obtained in our convention can be related 
to eigenvalues solved in standard 2x2 problem by shifting the phase of neutrinos and 
making the reduced hamiltonian traceless. 

We are ready to solve the evolution problem in ([T]). Note that H{x) can be re- written 

as 

H{x) = H + 5V{x), 5V{x) = V{x) - V. (6) 

We first solve the evolution by H and obtain the contribution of 5V using perturbation 
in 5V. Keeping result of first order in 5V we obtain 

i^{L) =M(L)V^(0), (7) 

M(L) =U^e-'^\l-iC) Ul, (8) 
where C is a 3 x 3 matrix accounting for the non-adiabatic transition: 

C= dx e'^''Ul5V{x)U^e-'^\ (9) 

It is clear that = C holds. One can see that 

Cn = t dx{Ul5V{x)Um)rj = ^. (10) 

C,k = ! dx e'^%Ul5V{x)U^),u. 3 ^ k. (11) 

Eq. ffTOl) is guaranteed by Eq. ([3]). \Cjk\ ^ l(j 7^ k) should be satisfied if this is a 
good perturbation theory. One of the virtues of this perturbation theory is that Eq. ffTOj) 
guarantees that the oscillation phase is correctly reproduced. 

In the next two sections we will discuss in detail that we indeed do expansion in 
small quantities when computing Cjk{j 7^ k) and correction of the second order (9(C^) 
is further suppressed by these small quantities. Hence we can convince ourselves that we 
are dealing with a perturbation theory. Before making detailed discussions on it we note 
that it is a perturbation partly because density changes mildly in the mantle or in the 
core of the Earth. In the mantle or in the core | |/| |^| | ^ 0.3. This helps in improving 
the perturbation theory. 
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When L > 10690 km neutrinos cross the core of the Earth. Large density jump between 
the core and the mantle can cause problem to the simplest version of the perturbation 
theory. To improve the approximation we use average potentials for parts of trajectory 
in the core and in the mantle separately. So the evolution matrix can be written as 

M = M3M2M1, (12) 

where M2 is the evolution matrix in the core and Mi^a are evolution matrices in the 
mantle. They are 

=[/™e-^^(^-^'-)(l-zC,) f/L, ^ = 1,2,3 (13) 

where Lq = and = L. x = Li is the point where neutrinos cross from the mantle 
to the core, a; = L2 is the point where neutrinos come out of the core to the mantle. Aj 
and Umi are the engenvalues and mixing matrix in the region x G [Lj_i,Lj]. They are 
computed using average potential 

V, = , \ / " dx Vix). (14) 

Ci is computed using Eq. in the region x G Because of the approximate 

symmetric density profile of the Earth we have 

Vs ^ V,. 

This strategy to treat evolution of core crossing neutrinos was also noticed in Ref. [T5] . 
In the present paper when we do computation using the perturbation theory we always 
use Eq. (1121) for core crossing trajectories which have L > 10690 km and use Eq. for 
trajectories crossing only mantle which have L < 10690 km. 
For anti-neutrinos the hamiltonian is 

H\x) = H^ -V^{x), (15) 

where V^{x) is the transpose of the potential term V{x) introduced in Eq. ([1]). is the 
transpose of Hq. The problem of neutrino evolution can be solved using average potential, 
similar to the program used for neutrinos. 

The probability of neutrino oscillation is introduced using the evolution matrix: 

Pii^k ^ i^i) = \Mikix = L)\\ P(z7fc ^ ui) = \Mikix = L)p, (16) 

where Mki is the evolution matrix of anti-neutrinos ip = diag{i/e, t'/x, J^r}: 

^P{L) = M(L)^(0). 
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Figure 1: P^u^ — > fe) versus energy for the case with standard matter effect. Left panel is for 
L = 3000 km and right panel is for L = 6000 km. Am^i = 8. x lO'^ eV^, Am^s = 3. x lO'^ 
eV^ sin^ 26123 = 1, tan^ (9i2 = 0.41, sin^ 26*13 = 0.01, ^13 = vr/G. PREM density profile is used 
for computation in this figure and all remaining figures in this paper. 

Using the oscillation probability we can define observables of CP violation and of time 
reversal asymmetry. For example we define 
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One can see in Eqs. ([1]) and (|T5l) that Acp does not contain pure information of fun- 
damental CP violating parameter. Because the medium contains only matter but not 
anti-matter Acp is not zero even if CP violating phase 613 in U is zero. It is clear 
that powerful electron neutrino beam is needed to observe the time reversal asymmetry. 
Intense electron neutrino beam is probably available using muon storage technology in 
future experiments. 

3 Oscillation with standard matter effect 

In this section we consider oscillation of neutrinos with standard matter effect. In this 
case we express the potential, mixing matrix, etc as 

V{x) = Vs{x), V = Vs, H = Hs, H = Hs, Um = Us, A = As. (19) 
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Figure 2: At versus energy for the case with standard matter effect. Left panel is for L = 3000 
km and right panel is for L = 6000 km. Parameters are the same as in Fig. [TJ 



The potential term of standard matter effect is 

= diag{V;(x),0,0}. (20) 

Ve = \/2GFNe is the potential in matter. A^e is the electron number density in matter. 
GiT' is the Fermi constant. Effect of standard NC interaction is universal for three types 
of neutrinos and is neglected for neutrino oscillation in Earth matter. Vs., Hg, Hg, Ug, Ag 
are defined using Eqs. Q, ([1]), ([2]), (jl]) and ([5]). Cjk is obtained as follows: 

Cjk = t dx e^T- {Ug):^{Ug)ek 5Ve{x), j ^ k, (21) 
Jo 

where 

dVeix) = Veix) - K, K = 7 / dx V,ix). (22) 

Jo 

For trajectories crossing mantle only the evolution matrix is 

Mg{L) =Uge-'^^{l-tC)Ul (23) 

For core crossing trajectories (L > 10690 km) the evolution matrix is generalized using 
Eq. ([12]). 

Mg{L) = MgsMg2Mgr. (24) 
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Ms2 and Msi,s2 are evolution matrices in the core and mantle separately. They are com- 
puted to the first order in 6V and using average potentials in the core and the mantle 
separately, as explained in section [21 

For comparison we also show the results computed using analytic formula: 

Mas{L) = Us e-^^^ Ul (25) 

where non-adiabatic contribution is not included. This formula is not generalized for core 
crossing trajectories as done in Eq. (12^ . 

In Fig. [T]we plot P{v^ —>■ Ve) versus energy. We show results of numerical computation, 
results computed using Eq. fl23|) and the result computed using analytic formula Eq. fl25l) . 
The lines labeled with "analytic with average Ve" are computed using Eq. fl25|) . The lines 
labeled with "first order in perturbation" are computed using Eq. ( 1231) . They are plots 
for L = 3000 km and 6000 km. One can see that results computed using the perturbation 
theory reproduce precisely the phase and the magnitude of the oscillation pattern. One 
can see that for these two baselines the analytic results computed using Eq. (1251) give a 
quite good approximation to the oscillation pattern. Actually one can not see difference 
in the left panel {L = 3000 km). In the right panel (L = 6000 km) one can see that 
oscillation phase is correctly reproduced by the analytic result but the magnitude around 
peaks is not precisely reproduced. There are some small differences. 

In the right panel height of the first peak of the right side is given by sin^ 26'^. One 
can see that it is much larger than the magnitude given by the vacuum value sin^ 2^13. 
This is because the position of the first peak is close to the region of 1 — 3 MSW resonance 
which has energy range ii^ ~ 7 — 10 GeV. Magnitudes of the second and third peaks are 
much smaller. This is because the second and third peaks are away from the resonance 
region and the vacuum value is dominant. In Fig. 121 we show plots for the time reversal 
asymmetry versus energy. Again in the left panel the analytic result using Eq. (1251) 
is precise. In the right panel there are some small differences between numerical and 
analytic results. Results computed using the perturbation theory, i.e. using Eq. (1231) . are 
in remarkable agreement with the numerical results. 

The analytic formula is a better approximation for shorter baseline. In Figs. and 
dl]) we plot Pii^fi — > i^e) versus L, the length of baseline. For L > 10690 km the lines 
labeled with "first order in perturbation" are computed using (1241) . Left panels in Figs. 
and dl]) are for E = 5 GeV and right panels are for ii^ = 15 GeV. One can see clearly 
that the analytic result is a very good approximation for L ^ 6000 km. It's no longer 
precise for L ^ 6000 km. In all these plots the results computed using the perturbation 
theory are always in remarkable agreement with that of numerical computation. 

In Fig. iniandiniwe also compare computations on P(z/^ — > I'r) and P(z/e Ur). We 
see that -P(f^ Vt) is very well approximated by the analytic result computed using Eq. 
(1251) except for L > 11500 km. For P{ue i't) there are some differences for L > 6000 
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Figure 3: ^(f^ — > fe) versus L the length of basehne for the case with standard matter effect. 
Left panel is for E = 5 GeV and right panel is for E = 15 GeV. Neutrino parameters are the 
same as in Fig. [TJ 




Figure 4: At versus L the length of baseline for the case with standard matter effect. Left 
panel is for E = 5 GeV and right panel is for E = 15 GeV. Neutrino parameters are the same 
as in Fig. [TJ 
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km. Differences of these magnitude are negligible in Piiy^ Vt)- Again we see that 
the results of first order in the perturbation theory are in remarkable agreement with the 
results of numerical computation. 

One can understand the success of this formulation of neutrino oscillation by noting 
that indeed we are expanding in small quantities and we are dealing with a perturbation 
theory. In the following we show that Cjk{i 7^ k) computed using Eq. (|2T|) is suppressed 
by some small quantities. Hence in (1 — iC) condition 

\C,k\ < 1 (jV k) (26) 

can be satisfied in the first order. The second order result of order 0{C'^) is further 
suppressed. Thus we can be confident on the perturbation theory. Note that another 
condition is that density changes mildly in the mantle or in the core, as already emphasized 
in section ([2]). (!26|) is explained as follows. We first consider the case with Am|^ > 0. 

i) For 0.5 GeV < E < 7 - 10 GeV, < K < holds and this is the region 
below 1—3 MSW resonance. In this range of energy the eigenvalues of H are approximately 

cos^ 6i2,Ve,^^) in the limit that correction of 0(sin6'i3) and 0{Aml^/{2EVe) ) 

are neglected. Small correction of 0{ ^^^ ) has been neglected in larger eigenvalues. Hence 
the first entries of H are changed to the second entries. We have 

Am^ 1 

sin2^^3 ^ sin 2^13, cos^™ ~ ^tI^tt sin 2^12. (27) 

4:E Ve 

So we get \{Us)e3\ = sin^™ < 1 and {Us)ei = cos ^^2 cos < 1. Cjk oc (t/,)^([/,)efc(j ^ 
k) is suppressed either by {Us)e3 or by {Us)ei- 

In this range of energy Cjk is suppressed by small quantities sin 6*13 and Am'^i/ {2EVe). 

ii) For E' > 7 — 10 GeV, Ve > holds and this is the region above 1 — 3 MSW reso- 
nance. In this range of energy the eigenvalues of H are approximately (-^^ cos^ 6^12, Ve) 
in the limit that corrections of 0(sin6'i3) and 0{Am\^/ {2EVe) ) are neglected. In this 
case the first entries of H are changed to the third entries. We can get (see appendix) 

Am^ 1 

cos^ra - ^TT^tF sm2^i3 « 1. (28) 

4-& Ve 

So we have {Us)ei = cos 6*5^ cos 6^^^ ^ 1 and {Us)e2 = cos 6'^ sin 6'^ <^ 1. Cjk is suppressed 
either by (t/Jei or by (f/s)e2- 

In this range of energy Cjk is suppressed by small quantity Am3^/(2£'V^) sin 2^13. 

iii) For ~ 7 — 10 GeV, condition Ve = cos 26*13 can be satisfied and the MSW 
resonance can happen. In this case cos 9^ is expressed using Eq. ( ETl) . C12 and C13 are 



9 




E (GeV) E (GeV) 



Figure 5: Plots for the case with standard matter effect. Left panel, ^(t'^ Vt) versus energy, 
L = 6000 km; right panel, Pive ^r) versus energy, L = 6000 km. Other parameters are given 
in Fig. [H 




Figure 6: Plots for the case with standard matter effect. Left panel, P{Vfji — > Vt) versus L, 
E = 5 GeV; right panel, P{i'e i^t) versus L, E = 5 GeV. Other parameters are given in Fig. 

m 
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suppressed by it. But C23 is no longer suppressed by sin^^ or cos 6*^ since 6*]^ ~ 7r/4 if 
resonance happens. 

One can see that C23 is suppressed by the resonance condition itself. In the resonance 
region the second and third mass eigenstates get almost degenerate: 

— (A^ - A^) ^ J sin2^i3 + (29) 

It is clear that in the resonance region 

A3- A2 

is a small number. So 

C23 = e-^T^ ["^ dx e-**(^) iUsrAUs)e3 
Jo 

= e-^T^/" dx iUs):,{Us)e3SV,ix)[l-t(Pix)-l(l)\x) + ■■■]. (30) 
Jo ^ 

The first term in the bracket gives zero after integration. The second term gives zero after 
taking into account the fact that 5V is approximately symmetric in the Earth. Since major 
contribution is from term it is not hard to see that IC23I is indeed suppressed by small 
numbers. 

In this range of energy Cjk is suppressed by tSTn^^ji^EV^ and Am\^L/{2E) sin 2^13. 

If Amg]^ < 0, eigenvalues of H are roughly (^^, cos^6'i2, K)- It is obtained by 
interchanging the first and third entries of H. Similar to the discussion on point ii) we can 
get cos ~ sin26'i3|Am3]^|/(2|Am|;^| + 4EVe). It is easy to see that we are expanding 
Cjk using small quantities {Us)lj{Us)ek{j 7^ k) which is suppressed by cos 6*^. 

For anti-neutrino the MSW resonance happens for Am|j^ < 0. Discussions on the 
perturbation theory closely follow the discussions for neutrinos above. Similarly, one can 
show that Cjk{j 7^ k) is expanded using small quantities. The second order of perturbation 
theory is of order (9(C^). Hence it is further suppressed. 

We summarize that we do expansion using quantities \{Us)lj{Us)ek\{j 7^ k) in the 
perturbation theory. Except in the MSW resonance region these quantities are small 
in energy range E ^ 0.5 GeV. In the resonance region the resonance condition itself 
guarantees the validness of the perturbative expansion. 

4 Oscillation with non-standard matter effect 

In this section we extend the discussion to case with non-standard matter effect. Previous 
works on non-standard matter effect include Refs. [2H ES] • In this case we express the 
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potential, mixing matrix, etc as 

V{x) = Vn{x), V = K, H = H^, H = H^, Ur,, = U„., A = A„. (31) 

Physics beyond the Standard Model can give non-standard four fermion interactions 
such as q'~ffj,q i^kl'^^^i and 67^6 u^'j'^i^i- Neutrino evolution in matter is modified by these 
terms. In general the potential can be written as follows 

/ Km KrX 

Vn{x) = diag{K, 0, 0} + V V , (32) 

\Vre Km VrrJ 

where Ve = \/2GFNe is the potential with standard charged current interaction, Vki is 
from non-standard NC interaction. VJ^ = Vm because the Hamiltonian is hermitian. x 
dependence in Vki has been suppressed in Eq. (!32|) . V^g has been made zero in our 
convention. This is achieved by shifting the phases of neutrinos: z/; — > e~*-^ '^^ ^""J^i, 
Vki Vki + Vee- In this convention Vki is 

Vki = \/2Gf [{fkl - fee)Ne + {Qh " ^ee)A^p + {hkl - K,)Nnl 

= Ve [{fkl - fee + Qkl - Qee) + {hkl - /lee) A^n/A^e] , (33) 

where fku Qki and hki are the dimensionless strengths of non-standard four Fermion in- 
teractions \/2 fkl Gp e7^e VkY^h ^ 9ki Gp p-f^p Vkl^i^i and a/2 hki Gp n-f^n Vkl^i^i- 
Np and are number densities of proton and neutron in matter. The second line in Eq. 
fl33|l holds because of the equality A^e = Np in neutral matter. 
We can re-write Vn as 

e/xT • (34) 

^re ^Tfj, ^TT / 

where eki = Vki/Ve- As can be seen in Eq. f l33l) . eki is not a constant if Nn/Ne is not 
a constant in matter. This happens when chemical composition changes in matter. In 
the present paper we are not going to consider the possibility that changes in the 
neutrino trajectory. We set constant in the analysis. This possibility is achieved when 
non-standard interactions are from interactions of neutrino with proton or with electron. 

Four fermion interactions of muon neutrinos are well constrained by direct test, e.g. by 
NuTeV experiment [22]. Since eki is a coherent combination of strengths fki, Qki and hku 
constraints on the strengths of these couplings from collider or fixed target experiments 
can not be directly translated to constraints on eki- However, one can have a rough 
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Figure 7: -P(t'^ fe) versus energy for the case with non-standard matter effect. Left panel 
is for L = 3000 km; right panel is for L = 6000 km. = 0.01 e-^'^/^o^ ^ q_o4 e-*'^/^, 
e^^ = 0.01 e"*'^/^'^. Other parameters are the same as in Fig. [1] 




Figure 8: At versus energy for the case with non-standard matter effect. Left panel is for 
L = 3000 km and right panel is for L = 6000 km. = 0.01 e'^'^/^o^ = o.04 e'^'^/^, 
e^^ = 0.01 e^*'^/20_ Other parameters are the same as in Fig. [TJ 



13 



constraint on the magnitude of en- The present experiments can reach precision of about 
one percent [23]. Hence, one can induce that |e^e|, Ic/^tI ^ 10^^. /ee,er,Tr, 5'ee,er,Tr, ^ee.er.rr 
are not well constrained because there are no powerful electron neutrino and tau neutrino 
beams. Hence in our convention egr, ^^j,fj. and err are not well constrained by direct 
tests. Other constraints on eki come from neutrino oscillation experiments. Previous 
studies [26l [27] show that |e^^|, {errl ^ 10~^ and {eerl ^ 10""^. 

The trajectory dependent average potential and average Hamiltonian are 

Ve Km 

H^ = Ho + K, K = I V,, V , (35) 

where Ve is given in Eq. (l22l) . Vfc; is similarly defined. Mixing matrix f/„ diagonalizes Hn-, 
as done in Eq. ([5]). A = A„ is obtained after diagonalization. 

Keeping the first order result in 5V the evolution matrix is solved as 

M„(x) =Une-''^%l-iC)Ul (36) 
C,k = r f/xe*T-[([/„)*^. (f/„)^, 5re(x)+ V (f/„);. (f/„)tfc 5Kt(x)], (37) 

where 5Vst = Vgt — Kt- Cjj = holds as explained in section [2J For core crossing 
trajectories {L > 10690 km) the evolution matrix is generalized using Eq. f[T^ as follows 

M„(L) = M„3 M„2 M„i. (38) 

where M„i „2 and M„3 are evolution matrices in the mantle and in the core. They are 
computed to the first order in SV and using average potentials in the core and in the 
mantle separately, as explained in section [2l 

We also show result computed using analytic formula 

M^niL) = Un e-*^^ Ul (39) 

This formula is not generalized for core crossing trajectories. 

In Fig. [7| we plot Piv^ — > Ve) versus energy for the case with non-standard matter 
effect. We show results of numerical computation, results computed using the pertur- 
bation theory and the results computed using the analytic formula Eq. (15^ . The left 
panel in the figure is for L = 3000 km and the right panel is for L = 6000 km. In the 



''^Intense electron neutrino or anti-neutrino source is available at low energy {E ^ 30 MeV) from 
stopped muon decay or from reactor. But they are not enough to make strong constraints. 
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Figure 9: ^(f^ — > i^e) versus L the length of basehne for the case with non-standard matter 
effect. Left panel is for E = 5 GeV and right panel is for E = 15 GeV. Neutrino parameters are 
the same as in Fig. [TJ 




Figure 10: At versus L the length of baseline for the case with non-standard matter effect. 
Left panel is for E = 5 GeV and right panel is for E = 15 GeV. Neutrino parameters are the 
same as in Fig. [TJ 
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right panel one can see clearly the effect of MSW resonance in the first peak from right 
hand side, sin^ 2^55 is enhanced to about 0.1, much larger than the vacuum value 0.01 
chosen in the analysis. The results computed in the perturbation theory are in perfect 
agreement with that of the numerical computations. The analytic result is not always in 
perfect agreement with that of the numerical computation. In the right panel one can see 
some differences of the two computations. These differences can also be seen clearly in 
the right panel of Fig. [8] where the time reversal asymmetry is plotted. 

In Fig. iniwe plot -P(t'^ —>■ Ve) versus L the length of the baseline. In Fig. [10] the time 
reversal asymmetry is plotted. One can see clearly the remarkable agreement between the 
computation using the perturbation theory and the numerical computation. The analytic 
computation is a good approximation for L ^ 6000 km. 

We can show that indeed we are expanding in small quantities and we are doing a 
perturbation theory. First, the second term in the bracket of f l37l) gives small contribution. 
This is because \Vu\lVe ^ 0.1. Hence |(5Vfci|L = (|5Vfci|/|VA;/|)| V^zl-f/ is a small number. 
Second, the discussion on the first term closely follows the discussions in the previous 
section. In lower energy region ( 0.5 GeV E < 7 — 10 GeV) the discussion in the 
previous section is still valid since Vki does not change the mixing matrix very much. For 
energy > 10 GeV, V^i can be important. When Vki is important enough to determine 
mixing matrix we have (see appendix) 

cos ^^3 ~ \l\Vre? + Ke\VVe < 1- (40) 

Contribution of the first term is suppressed by {Un)ei,e2 which are proportional to cos^]^. 

We summarize that in the presence of non-standard matter effect the perturbation 
theory is valid because it is expanded using small quantities {Us)lj{Us)ek{,i 7^ k) and dVki- 

5 Discussions and conclusions 

In this paper we propose a perturbation theory which simplifies the problem of neutrino 
oscillation in the Earth. We perform analysis with three flavor of neutrinos. The pertur- 
bation theory is developed using the trajectory dependent average potential. The average 
potential is averaged along the neutrino trajectory in the Earth. So it depends on the 
trajectory of neutrino. The effect of non-constant density profile in the Earth, i.e. the 
non-adiabatic contribution, is carefully included in the first order of the perturbation the- 
ory. This perturbation theory is generalized using Eq. (fT2l) for core crossing trajectories, 
i.e. for L < 10690 km. The problem of neutrino oscillation is substantially simplified 
using the perturbation theory presented in the present paper. 

Using the perturbation theory we study neutrino oscillation in the Earth for cases with 
standard matter effect and with non-standard matter effect. It is shown that for both cases 
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the formulation presented gives a precise description of the neutrino oscillation. We study 
observables of flavor conversion and of time reversal asymmetry in neutrino oscillation. 
We find remarkable agreement between result of computation using our perturbation 
theory and that of numerical computation. This is a nontrivial agreement. To see that 
it is nontrivial one can try another perturbation theory. For the case with non-standard 
matter effect one may rewrite the hamiltonian as follows 



One may first solve evolution problem of Hi and treat H2 as perturbation. The first order 
result of this perturbation theory is a good approximation for < 5 GeV. Apparently it 
will not be good for higher energy neutrinos when non-standard matter effect becomes 
important. It is not able to reproduce the time reversal asymmetry in high energy region 
which is dominated by non-standard matter effect. To improve the perturbation theory 
one may go to higher order and make the description more complicated. Apparently this 
perturbation theory is not a good choice compared to one given in the present paper. 

We also show the results of analytic computation, i.e. the computation using average 
potential but without non-adiabatic contribution. It is shown that the analytic result 
gives a very good approximation to the oscillation for L < 6000 km. It means that the 
problem of neutrino oscillation can be greatly simplified for L ^ 6000 km and E ^ 0.5 
GeV. For this range of parameter space matter effect in a fixed baseline is described by 
a constant potential term. For standard matter effect, in particular, it says that matter 
effect is a one parameter fit. 

The perturbation theory given in the present paper is valid for ^ 0.5 GeV. Previous 
works [31 S] cover the energy range E ^30 MeV. They are all energy ranges away from 
the 1 — 2 MSW resonance. It is a question to find a compact and simple formulation of 
neutrino oscillation in the Earth in the energy range 30 MeV ^ E ^ 500 MeV for which 
the 1 — 2 MSW resonance can happen in the Earth. 

We perform analysis using the PREM density profile [TB]. It is a symmetric density 
profile. We also assume eki is constant in considering non-standard matter effect. These 
assumptions are not generally true. For more general density profile further research 
is needed to understand whether the perturbation theory is as good as analyzed in the 
present paper. We expect that this perturbation theory is still a good approximation in 
more general cases. Research on neutrino evolution with more general density profile will 
be presented in further publications. 
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H = Hi + H2, Hi = Ho + dmg{Ve,0,0}, 




) 
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Appendix 

In this appendix we show that in large Ve hmit cos 6^ is smalL We write 

if = if" + diag{K,0,0}, 

where 

= Ho, 

for the case with standard matter effect, and 

/ Ve^ Ver\ 
H^ = Ho+\ V^e V^r , 

\Vre Km K./ 

for the case with non-standard matter effect. 
We consider the eigenvalue problem of H: 

HX = XX, X = {xi,X2,X3f. 

In the large Ve limit one of the engenvalues of H is roughly Vg. So at zeroth order A 
We solve the eigenvalue problem using perturbation in H^/Vg-. 

A = Ao + Ai + As + ■ ■ ■ , X = Xo + Xi + X2 + • • • 

So to first order in H^/Ve we get 

diag{V;, 0, 0} Xo = Ao Xo, 

Xo + diag{K, 0, 0} Xi = Ao Xi + AiXq. 

We get from Eq. f H6i) 

Ao = K, Xo = (l,0,Of. 
Putting Eq. (gH]) into Eq. (gTD we get 

Ai = H^e, Xi = (0, H^jVe, KJVef. 
Writing X = (sin6'55e~*''i3, cos6'55 sin6'^, cos6'5^ cos 6'^)'^, we can get 

cos''e^,^{\Hlf + \H:f)lV^. 
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Using Eq. in the case with standard matter effect we can get 



cosOTs = ^^Vl-|(f/o)e3|^ \{UoU = ^-^szn29,,. (51) 

In the case with non-standard matter effect there are two possibihties. One possibihty is 
that effect of Hq is larger than the non-standard matter effect. Hence, Eq. flSTj) holds and 
we have small cos^^ in large Ve limit. The other possibility is that non-standard matter 
effect is large enough to determine the mixing matrix. So we have 

cos ^^3 = ^JKe?+\Vre\VV,. (52) 

We conclude that cos^^ is a small quantity in large E limit for both the case with 
standard matter effect and the case with non-standard matter effect. 
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